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Abstract 

Improved computation of the dielectric function considering excitonic effects and long wave¬ 
length is performed and compared with the nearly free electron band approximation, similarly 
with the Penn’s model case. New expressions for the real and imaginary part of the dielectric 
function are presented and the real part compared with the Penn’s result. The obtained functions 
satisfy the Kramers-Kronig relations, in contrast with earlier results in the literature. In addition, 
our improved dielectric function presents a coefficient of 2/3 for small gap approximation (different 
from the value of 1 in the original Penn model) is very close to the value 0.62 obtained in [Can. 
J. Phys.53,(1975) p.2549] from pure numerical procedures. The obtained dielectric function also 
is used in a rough and stimative analysis of the metal-insulator transition in molecular hydrogen 
being the critical densities determinated near the experimental values for the hydrogen coming from 
other approach. The approximated expressions and critical values are given and the usefulness of 
the rough methods involved in the determination of the critical points briefly discussed. 
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I. THE COMPLEX DIELECTRIC FUNCTION 


As it is well known, the refractive index and the energy gap of semiconductors represent 
two fundamental aspects that characterize their optical and electronic properties. Penn in 
1962 [1] proposed a simple model for an anisotropic semiconductor with electrons in a sphere 
of momentum space and by an isotropic energy gap. Penn’s model, despite its simplicity, was 
and is the basic ingredient of several physical applications to diverse materials and different 
problems of great interests as the insulator-metal transition, exciton condensation and due to 
the isotropy of the system the model can be successfully applied to a liquid or an amorphous 
semiconductor. In this short report we improve the computation of the dielectric function 
going beyond the simple semiconductor considering the possibility of formation of excitons 
into the system. The obtained dielectric function (that fulfils the Kramers-Kronig relations 
[4]) is briefly analyzed and compared with the original Penn’s one in its real (static) part in 
order to see the main differences and how the excitons modify the analytical properties of 
the dielectric function (critical points, singularities, zeroes, etc.) . 

Our starting point is the expression for the complex dielectric function in the reduced 
zone [2] 
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where fo ( E ki) is the distribution function for the reduced wave k and the band l, Q :volume 
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of the solid, | k,t) is the Bloch wave function that satisfy H 0 | k,t) = E^ t \k,t) where H 0 is 
the Hamiltonian for an electron in the unperturbed lattice, e (ca, q) in the extended zone 
scheme takes the following form: 
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where the distribution functions were changed by the respective zone occupation numbers. 
The complex dielectric function is defined, however, as e (to, q) = (ca, q) — ie 2 (to, q) 
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We are interested in t\ (0, q) and e 2 (to, q ) in the Lim q —* 0 and to obtain closed expressions 
for these quantities from the above expressions (similarly as in the case of the Penn model 
for semiconductors). 

The model is equivalent to the NFE (nearly free electron) model isotropically extended 
to 3 dimensions. The energies of the 2 bands are given by 
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k = k — 2 kpk, and 2 kpk is the Fermi wave vector. Let us consider 


the wave functions (r) = ^ 


0 ik-r 

72 


H k _ 

1/2 


L Mf) 


where = 


± _ ( Eg-Ex ) _ A 


2 ET 1 = ifar, is the volume 

fc k 


and the signs — (+) corresponds to the first (second) Brillouin zone. We suppose now ( 
following the Srinivasan condition [3]) that k and k + q have the same reciprocal vector 
—2 kpk = G (e.g. Umklapp process). Substituting into the general expressions (3,4), they 
depend now only on k,q 
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and 
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Formally, we need to consider the following assumptions due to the specific physical problem 
that we are involved in: 

i) u is restricted to positive frequencies, then, N-% — 1 iff k G l st Brillouin zone, N-% = 0 iff 
k G elsewhere , 

ii) to perform the integrals, as usual, we assume a real crystal. 

The evaluation of the matrix elements yields: 
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when q —> 0 the leading order of the matrix elements yields 
\(k\e- iq - r \k + q)\ 2 = (k\e~ iq - r \t + q^ ' « 


fc+<?/ 

( 8 ) 


(9) 


cos 2 6 (A/ AE F y (q/k F ) z _ p 2 x 2 £ 2 


2 ’ 


4 [(1 - £;/£; F ) 2 + (A/4£ F ) 2 ] 4 C ^ 2 + P 2 ) 

in the standard notation: 9 is the angle between k and g, Ep = , A = E g — E x and 

A/4 Ep = p. If the normal process does not contribute to €\ (0) or 62 (E), from previous 
expressions we have 
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where the limit q —> 0 has been taken. The sum over k (l s< Brillouin zone) is converted in 
an integral 
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and 
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A -p(E 2 - A 2 ) 


1/2' 


E = hut, 


2 E (E 2 — A 2 ) 1/2 

the second equation holds for A < E < Ep (1 + p 2 ) 1 ^ 2 and e 2 = 0 otherwise. 
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For Penn’s case 
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clearly in this case there is a double pole (singularity)when E x —* E g and we see the evident 
mo dif ication to the shape and critical points of the dielectric function due to the presence 
of excitons in the system. A full analysis considering jointly the quantum held theoretical 
part of the history will be present elsewhere [4], For the imaginary part, when A 2 —>■ E 2 , 
e 2 (E) —> oo, and when A = 0, e 2 (E) = 0. In the next Section we will see the consistency of 
the above results making use of the Kramers-Kronig relations. 

II. KRAMERS-KRONIG RELATIONS 


Here we will corroborate the results of the previous paragraph by mean of the use of the 
Kramers-Kronig relations (following the same approximations as in previous Section). To 
begin with, from the relations between eiand e 2 we have 
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In the above equation the principal value symbol can be dropped: the delta function in 
e 2 ( E') restricts its non-null range to A < E < Ep (1 + p 2 ) ~. The crucial point now is the 
evaluation of the delta function in e 2 up to the second order in (/.From eq. (5) (previously 
to make the limit q —> 0) the energy difference is 
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with the definitions that were introduced before. 

Now we convert the sum in k in the first Brillouin zone in an integral obtaining 
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where the lower limit in the p integration is zero. We must keep in mind here that q is small 
but not zero, then, we must retain the relevant order in q. We write the delta function as 
8 (f 0 * 0 ) = twj8 (.x - x 0 ) with 
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As in several references based in the analytical models of semiconductors[ 6 ], we investigate 
the case E = A 

i) At rj = 0 , / (rj) > 0 (because the limit q —> 0 has not been taken.) 

ii) For rj ^ 0, 4 Ep (p 2 + p 2 ) 1 ^ 2 — A is positive and always dominates for sufficient small 
p (or q). Therefore, / (p) > 0 for E — A and the delta function cannot be satisfied. 

This clearly means that the point E = A has no contribution to 62 (i?). Putting / (p) = 0 
we determine p 0 in the limit g —» 0 in a similar manner that earlier works 
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The expression for p 0 is valid for all E>A + O (£ n ), 0 < n < 2. The region of invalidity 

vanishes in the limit £ —» 0 except in the precise point E g — E x . Finally we have 
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putting all these ingredients in the expression for 62 (-F) we obtain 
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we immediately obtain 
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arriving to the same result as in the formula (13) of the previous Section. In an analogous 
manner and straighforwardly by mean eq. (16) and (24) we arrive to expression (12) as is 
required by consistency. 


III. DIELECTRIC SCREENING AND INSULATOR METAL TRANSITION 

As it is well known, one way to study the transition coming from the insulating phase, 
is supplied by the static dielectric function e (0) of the molecular crystal which will be 
singular at the transition point because, as is expected, the transition to the metallic state 
is preceded with softning of some characteristic low lying excitation. At normal pressures 
and low temperature the ’’Frenkel exciton” is formed for a tightly bound electron hole pair. 

In sharp contrast, when the pressure is considerably augmented to some critical value 
on the molecular crystal, the intermolecular interaction increases and the electron and the 
hole move further away from each other. Then, we can expect that at such high density in 
the non-conducting phase the excited state is a (much extended) Wannier exciton. In this 
case, the detailed lattice structure of the crystal is not very important to the pair and the 
common procedure is to simulate the lattice effects introducing effective electron and hole 
masses and screened Coulomb interaction ~ 1 /er e h to describe an isolated Wannier exciton. 
Then, the ground state of the energy of the electron hole pair is simply -E x (centre of mass 
kinetic energy term avoided) where 


Ex — 2 [Ry\i ( 25 ) 

/ / i C 

is the Wannier exciton binding energy. However, // is the reduced mass, m is the bare elec¬ 
tron mass and e is the dielectric static function of the medium. In some earlier attempts, 
was claimed that is possible to use roughly this information to construct the correspond¬ 
ing effective model of the static dielectric function e\ (0) of the medium on the low lying 
excitations in the case of compressed molecular hydrogen. We will follow such attempt in 
order to test it, due that our expression for the static dielectric function has a more detailed 
structure than the Penn’s case (that is also the case of ref. [10]). 
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As we have computed, the static dielectric function (15) near the metal insulator transi¬ 


tion Eg ~ E x becomes up to the lowest order 
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that clearly differs from the Penn case (and for the case of the ref. [10] for E g E x by a 2/3 
factor into the second term. In practice this equation needs to be solved self-consistenly: it 
is easy to solve the equation (26) if we know how /i and E g simoultaneusly depend on t\ (0) 
provided that the density limit of the existing electron-hole pair is maintained. Although 
the difficulty to find a full self-consistent solution for molecular hydrogen of 6\ (0) certainly 
exists, it is possible to obtain some estimates of the electric function near the insulator-metal 
transition. Consequently, our argument is that the exciton radius X x cannot be higher than 
a certain characteristic length (g* ) 1 . Here q FF is the critical Thomas-Fermi screening 
constant associated with a screened Heitler-London H 2 molecule immerse in an electron gas 
at fixed mean density prior dissociation: if the mean electron density of the gas is such that 
q TF > q* F the screened H 2 molecule turns to be unstable. Using variational, perturbative 
and other methods [9] diverse values for g* F (see Table I below) were found. These values 
can be used as upper limit for that is in the transition region 
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If the neighborhood of the insulator-metal transition E x varies slowly as function of the 
density in comparison with E g we can substitute the result (27) back into (25) to obtain 


E ~ E* 



(28) 


We see that near the insulator-metal transition E* is saturated at this critical value. Using 
Friendly-Ashcroft band structure calculation [7] of the E g for the molecular hydrogen under 
pressure together with the combination of (28) and (26), we get the following values of E* 
for the diverse methods [9] used to obtain q TF 




METHOD q TF E*[Ry\ r* 

Variational molecular orbital 1.197 0.067 1.56375 

electron pair 1.166 0.068 1.565 

Perturbative 1.000 0.080 1.58 

Other(average ) 1.190 0.0672 1.564 

Notice the lack of precision in r* despite the different values of q TF , E*[Ry\ obtained by 
the different methods listed into the table above: evidently attempts as these [ 10 ] must be 
accompanied with other ones, looking for the microscopic/topological side of the exciton 
systems. 

IV. CONCLUDING REMARKS 

In this paper an improved derivation of the dielectric function based in a similar analytical 
model as the Penn proposal was presented, but now considering excitons in the system. 
Also in our calculation the Umklapp process as in [3] and long wavelength (q —> 0) was 
considered. The differences with Penn’s model in the case of the real part of the dielectric 
function was given showing how the presence of excitons in the system modifies its behaviour. 
These results, being of a theoretical analytical character, were corroborated by mean of the 
Kramers-Kronig relations arriving independently to the same expressions for 62 ( E) and 
61 ( E ) up to the required order. Also it is important to note that our pure theoretical result 
presents a coefficient of 2/3 for small gap approximation (different from the value of 1 in the 
original Penn model) is very close to the value 0.62 numerically obtained by R. D. GRIMES 
and E. R. COWLEY in [11], Using the improved dielectric function and making some 
standard assumptions about the behaviour of the exciton binding energies near a saturation 
value before the dissociation, we showed that e\ (0) diverges at the transition point. This fact 
is in complete agreement with the experimental result of [ 8 ] and it stresses a discontinuous 
nature of the metal-insulator transition in solid hydrogen at low temperature under pressure. 
However same computation for different values of q TF , E*[Ry] obtained by the different 
methods (listed in the Table) present a lack of precision showing the roughtness of this type of 
approaches. However, the rough analysis of such a type of transitions must be complemented 
from a microscopic point of view. Consequently, it will be of particular interest to use in 
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near future the model presented here together with a ” Keldysh-like” approach as [5] in order 
to analyze the metal insulator transition. 
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